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Abstract. We construct Green's functions for divergence form, second order parabolic 
systems in non-smooth time-varying domains whose boundaries are locally represented as 
graph of functions that are Lipschitz continuous in the spatial variables and 1 /2-H61der 
continuous in the time variable, under the assumption that weak solutions of the system 
satisfy an interior Holder continuity estimate. We also derive global pointwise estimates for 
Green's function in such time-varying domains under the assumption that weak solutions 
of the system vanishing on a portion of the boundary satisfy a certain local boundedness 
estimate and a local Holder continuity estimate. In particular, our results apply to complex 
perturbations of a single real equation. 



1. Introduction 

Green's functions play an important role in the solution of elliptic and parabolic partial 
differential equations. There is a large literature on Green's functions of uniformly ellip- 
tic and parabolic equations in divergence form. Green's functions of elliptic equations of 
divergence form with L M coefficients have been extensively studied by Littman et al. 11231 
and Gruter and Widman lfl2l ; see also Q OH HH . Recently, Hofmann and Kim lfl4l gave 
a unified approach in studying Green's functions for both scalar equations and systems 
of elliptic type; see also (8). For parabolic equations, Aronson [T) established two-sided 
Gaussian bounds for the fundamental solutions of parabolic equations in divergence form 
with Loo coefficients; see also 0[6][9][T3][T8l|28] an d references therein for related results. 
In a recent paper by Cho and the authors (U , we proved that if weak solutions of a given 
parabolic system satisfy an interior Holder continuity estimate, then the Green's function 
of the system exists in any cylindrical domain. In the scalar case, such an interior Holder 
continuity estimate is a consequence of Nash 1251 and Moser l24l . and also such an esti- 
mate is available for weak solutions of a system if, for example, its principal coefficients 
are VMO in the spatial variables. However, the construction of Green's function in JU 
heavily relied on the results by Ladyzhenskaya and Ural'tseva that are available only for 
cylindrical domains. In another recent article by Cho and the authors 0, we demonstrated 
how one can derive global pointwise estimates for the Green's function in a cylindrical 
domain by using a local boundedness estimate and a local Holder estimate for the weak 
solutions of the parabolic system vanishing on a portion of the boundary. 

The aim of this article is to give results similar to those of |U[5) for a class non-smooth 
time-varying domains whose boundaries are given locally as graph of functions that are 
Lipschitz continuous in the spatial variables and 1 /2-H61der continuous in the time vari- 
able, which hereafter shall be referred to as time-varying H\ domains. There are many 
papers dealing with parabolic equations in this type of time-varying domains. Lewis and 
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Murray 11211 considered a domain in R x R of the form {(t,x): xr > f(t, x')}, where 
x' = (x , ■ ■ ■ ,x ct ~ l ) and / is a function that is Lipschitz in x' and whose f-derivative of 
order 1/2 belongs to BMO, which is slightly stronger than f € Hi, The non-cylindrical 
domains considered by Hofmann and Lewis in their important paper ||15| on L2 boundary 
value problems for the heat equation are also included in time-varying Hi domains; see 
also lfl6l l27l l29l . Brown et al. [3) investigated weak solutions of parabolic equations in 
time-varying Hi domains and proved the unique solvability of Dirichlet boundary value 
problems. We shall in fact utilize their result in constructing the Green's function in time- 
varying Hi domains. In contrast, there is little literature on Green's functions of parabolic 
equations in non-cylindrical domains and to the best of our knowledge, there is no literature 
dealing with Green's function for parabolic systems of second order with Loo coefficients in 
time-varying Hi domains; see the remarks made in the last paragraph of the introduction. 

We denote by u — (u , . . . , u"') T a vector-valued function of d + 1 independent variables 
(f, x l , . . . , x d ) = (f, x) — X. We consider parabolic systems of second-order 

(1.1) Jz?h := /// - D a {A aP D p u\ 

where the usual summation conventions are assumed and A" 13 = A"^(X), for a, ft = 
l,...,d, are m x m matrices whose entries are L M functions satisfying the strong ellip- 
ticity condition; see Section F2~2l for the details. We emphasize that the coefficients are not 
assumed to be time independent or symmetric. We will later impose some further assump- 
tions on the operator Jz? but not explicitly on its coefficients. By a Green's function for the 
system Jl.lt in a time-varying Hi domain Q we mean an m x m matrix valued function 
G(X, Y) = G(t, x, s, y) which satisfies the following for all Y e Q: 

JzfGO, Y) = 6 Y (-)I„, in O, 
G(-, Y) = on @>a, 

where dy(-) is a Dirac delta function, I m is m x m identity matrix, and .^Q. denotes the 
parabolic boundary of f2; see Section 12.61 for more precise definition. In this article, we 
prove that if weak solutions of d 1 - 1 b satisfy an interior Holder continuity estimate, then 
there exists a unique Green's function in Q. and it satisfies some natural growth properties; 
see Theorem 13.11 below. Moreover, we show that the Green's function also satisfies the 
following familiar property: 

lim G(t, x, s, •) = 6 x (-)I m on oj(s) = {y e R d : (s,y) e Q). 

We also derive the following global Gaussian estimate for the Green's function in a time- 
varying Hi domain by using a local boundedness estimate for the weak solutions of ( II . lb 
vanishing on a portion of the boundary: For any T > 0, there exists N > such that for all 
X = (t, x) and Y = (s,y) in Q satisfying < t - s < T, we have 

N ( k\x - vl 2 

(1.2) \G{t,x,s,y)\< exp \-±^L 

(t - s) ' [ t—s 

where k > is a constant independent of T; see Theorem 13.1 II and Remark 13.131 In 
particular, the above estimate ( 11.2b holds in the scalar case (i.e., when m = 1) and also in 
the case of L x -perturbation of diagonal systems; see Corollary 14. H and Section |4~2l below. 
In fact, in such cases, a stronger estimate is available near the boundary. For any T > 0, 
there exists N > such that for all X = (t, x) and Y — (s,y) in O satisfying < t — s < T, 
we have 

(1.3) |G(^,v)|^(lA^niA^f * exp/ ^ 



\X-Y\ ?? ) \ \X-YW) (f-sW 2 



GREEN'S FUNCTIONS IN TIME- VARYING DOMAINS 



3 



where k > and jd e (0, 1] are constants independent of T, and we used the notation 
a A b = min(a, b\ \X - Y\g» = max( ^/\t=7\, \x - y\), and d(X) = inf{|Z - X\ Z e 3D,}. 
We show how to derive a global estimate like ( 11.31 for the Green's function in a time- 
varying Hi domain by using a local Holder continuity estimate for the weak solutions of 
(11. Il l vanishing on a portion of the boundary; see Theorem 13. 161 and Remark [3.191 As 
mentioned above, the estimate ( 11.3b particularly holds in the case of Loo-perturbation of 
diagonal system as well as in the scalar case; see Corollary 14. 4l and Section l4~2l 

The organization of the paper is as follows. In Section [2] we introduce some notation 
and definitions including the precise definitions of time-varying H\ domains and Green's 
functions of the system (II . lb in such domains. In Section [3] we state our main theorems 
and give a few remarks concerning extensions of them. In Section |4] we present some 
applications of our main results including applications to the scalar case, Leo-perturbation 
of diagonal systems, and systems with VMO, v coefficients. We provide proofs of our main 
theorems in Section|5]and some technical lemmas are proved in the appendix. 

Finally, several remarks are in order. In the scalar case, there are a few papers discussing 
Green's functions in non-cylindrical domains. However, we believe that even in the scalar 
case, our results give still new perspectives on Green's functions. In 11261 . Nystrom con- 
structed Green's functions in bounded time-varying H\ domains utilizing the fundamental 
solutions and the caloric measures, and in doing so, he made a qualitative assumption that 
the coefficients are smooth in order to have well-defined concept of solutions; i.e. to as- 
sume that all solutions are classical ones. The main drawback of this kind of approach is 
that it is not well suited to handle unbounded domains, especially domains with unbounded 
cross-sections such as the graph domains considered by Hofmann and Lewis fl5l . The 
novelty of our paper lies in presenting a powerful unifying method that establishes the ex- 
istence and various estimates for the Green's function of parabolic equations and systems 
with Loo coefficients in time- varying H\ domains including the graph domains. Also, even 
though we impose some conditions on the operator _Sf in the vectorial case, we do not 
make any smoothness assumption on its coefficients in order to assume that the solutions 
of the system are classical. Moreover, the treatment of Loo-perturbation of diagonal sys- 
tems is a unique feature of our paper and we believe that it could find some interesting 
applications in the complex perturbation theory for the Dirichlet problem of second order 
parabolic equations in time-varying domains. 

2. Notation and Definitions 

2.1. Basic notation. We mostly follow notation employed in Ladyzhenskaya et al. l20l . 
supplemented by that used in Lieberman l22l . First we use X = (f, x) = (t,x\ . . ., x' 1 ) to 
denote a point in R d+1 with d > I and we denote X' = (f, x') = (t,x 1 ,.. .,x d ~ i ) e W 1 so that 
X = {X'y). We also write Y = (s,y) = {s,y',y d ) = (!",/). We denote 

a A b — min(a, b), a V b — max(a, b) for a, b € [— oo, oo]. 

We define the parabolic distance in R d+1 and R d , respectively, by 

\X-YW = VF^Vlx-yl, \X'-Y'\& = y/\T^\V\ X ' -y'\, 

where | ■ | denotes the usual Euclidean norm, and write \X\gg = \X - 0\gs. We define the 
parabolic Holder norm as follows: 

\u\ M /2#-xi = M///2, w e + \u\o-n := sup + sup \u(X)\, (i e (0, 1]. 

XJeCI I A - / I at Xefl 
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By C M ' 2 ' M (Q.) we denote the set of all bounded measurable functions u on Q. for which 
\u\)i/2,ii;n is finite. We write D,u = D x iu = du/dx' and u, = du/dt. We also write Du = D x u 
for the vector (D\u, . . . , Dju). We use the following notation for basic cylinders in R d+l : 

Q(X ,R) = {X e R^ 1 : \X - X \&> < R}, 
Q-(X ,R) = {X = (t,x) eR d+1 : \X - X \& < R, t < t ], 
Q + (X ,R) = {X = (t,x)eM. d+l : \X -X \& < R, t > t }. 

We also use the ball B(xq, r) = {x e R d : \x- xo\ < r). For convenience, the parameter Xq 
(or xo) in the notation above is omitted if Xq — (or xq = 0, respectively). We use Q. to 
denote a domain (open connected set) in R d+l . For a fixed number fo, we write a)(to) for 
the set of all points (to, x) in Q, and write /(Q) for the set of all t such that co{f) is not empty. 
For -oo < to < t\ < oo, we denote 

Q(f ,/i) = {X = (t,x) € Q: t o < t < fi}. 

The parabolic boundary is defined to be the set of all points Xq e <9£2 such that for any 
s > 0, the cylinder Q-(Xq, s) contains points not in O. We define BO to be the set of all 
points X e such that there is a positive R with Q + (X ,R) c O and SO = \ BD.. 
We define the "time-reversed" parabolic boundary SPQ. to be the set of all points Xq e dCl 
such that for any e > 0, the cylinder Q+(Xq, s) contains points not in O. We also define 

Q[x, R] = Q n g(x, /?), ^Q[i, /?] = ^>Q n <2(x, /?), ^h[i, = n g(x, /?), 

and similarly £2 ± [X, R], £PQ. ± [X,R], and ^ 3 Q. ± [X, R]. Finally, we define distance functions 

d n (X) = d(X) = inffJF -X\&>: Ye d£l}, 
d n (X) = d~(X) = inf{\Y -X\&: Ye @>£l, s < t), 
d^(X) = d + (X) = M{\Y - X\g> : Y e s > t). 



2.2. Strongly parabolic systems. Let the operator Jzf be defined as in dl.ll i. We assume 
that the coefficient of «£? are defined in the whole space R d+l in a measurable way and that 
the principal coefficients A"^ with the components A a ? satisfy the strong ellipticity 

m d m d 

(2.1) 2 Z^ (X) ^^ v ZZlrf =:v lrf' VX€R" +1 , 

i,j=l a,fi=l i=l or=l 

and the uniform boundedness condition 



(2.2) 




VIef 



y=i 0-^8=1 



for some constant y e (0, 1]. The adjoint operator r _S? is defined by 

'Sfu = —Ui - D a {A ap Dpu), 

where A a/3 = (A 13 ") 1 '; i.e., A? 8 = A^". Notice that the coefficients A?? satisfy the conditions 
( 12.lt and ( 12. 2t with the same constant v. 
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2.3. Time-varying H\ domain. We shall say that O is a time-varying H\ graph domain 
in K.^ 1 if there is a function / = f{X') = f(t, x') from W 1 to R satisfying 

(2.3) \f(X')-f(Y')\<M\X'-Y'\^, Vr,F'6R d , 

for some constant M > so that Q is represented by 

Q = {X = (X', 6 R^ 1 : jc* > /(X')). 

We shall say that Q is a time-varying Hi domain in R d+l if 

i) 7(Q) = R and u>(t) is a bounded domain in W 1 for all f 6 R. 

ii) There are constants M and R a > such that for each e dQ., there is a function 
/ = f(X') = f(t, x') satisfying (12.3b . for which (after a suitable rotation of x-axes) 

Q n Q(X ,R a ) = {X e Q(X ,R a ): x d > f(X% 

2.4. Function spaces. For q > 1, we let L q {D) denote the classical Banach space consist- 
ing of measurable functions on O that are g-integrable. The space W2' (Q) denotes the set 
of functions u e L q {Q) with its weak derivative Du e L 9 (Q) having a finite norm 

We denote by W,' (fl) the Hilbert space with the inner product 

{u,v) w i,u a) := I uv+y I D a uD a v + I u,v,. 
1 Jn ^-f Jn Jn 

We define V2(^) as the set of all u e W^'^Q) having a finite norm |M|y,(n) defined by 



||w||y 2 (Q) := J \Du\ 2 dX + ess sup I u 2 dx. 

2 Jn tei(o.) Joj(t) 



The space Vj'^Q) is obtained by completing the set W^'^O) in the norm of V2(£2). Let 
E c Q and u be a V° ' (Q) function. We say that u vanishes (or write u — 0) on E if u is a 
limit in V^''(Q) of a sequence of functions in C™(Q. \ 2"). We define V°'\n) to be the set 
of all functions u in ^''(Q) that vanishes on SCI. 

2.5. Weak Solutions. For f,g a e Li j / oc (Q) m (a = 1, ...,d), we say that a is a weak 
solution of «5ftt = / + in Q. if m € V2(Q) m and satisfies 

(2.4) - f f AfB^Dj= f /y- f ijDafi, V0eC c °°(Q) m . 

Jn Jn Jn Jn 

We say that h is a weak solution of = / + Z) a g tr in Q if h e V r 2(Q)" i and satisfies 

(2.5) f h'>; + f A ? Dpu j Daft = f />' " f SaDaft, e C c °°(n) m . 
Jn Jn Jn Jn 

For if/ = if/ (x) € L\ i oc {u){ti))) m , we say that m is a weak solution of the problem 

Jz?h = / + D a g a in Q(fo,fi), u — on SQ(fo, fi), u = iff on w(fo) 
if m e ^''(Ofo, £i)) and satisfies for all t e /(£2(*o, fi)) the identity 

I u'ftdx- | u i (f>\dX+ \ AfD /3 u j D a (p i dX= \ fftdX 

Joj(t) Jn(r ,r) Jn(f ,T) U Jn(«o,T) 

- f g i a D a <p i dX+ f ifr^dx, V</>eC™(Q.(to,h)\SCl(t , 
Jn(/ ,r) Jaw 



fi)) m - 
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Similarly, for if/ = iff (x) e Lij oc (a>(ti)) m , we say that u is a weak solution of the problem 

'Jz?h = / + D a g a in Q(/o,/i), (i=0 on SQ(to,t\), u = ifr on u){t\) 
if u e V^fifo, fi)) and satisfies for all t e 7(fi(? , fi)) the identity 

I u i <f> i dx+ I u i 4> i t dX+ ( AfD /3 u j D a cp i dX= \ f'fidX 

- f g i a D a cf> i dX+ f tf/dx, V^eC c °°(0(^)\50(fo,/i)) m . 

2.6. Green's function. Let fi be a time-varying //[ (graph) domain in R d+1 . We say that 
anmXm matrix valued function G(X, F) = G(f , x, s, y), with entries Gjj(X, F) defined on 
the set {(X, Y) e fixfi : X + F}, is a Green's function of Jz? in if it satisfies the following 
properties: 

i) G(-, Y) e W°'/ oc (fi) and Jz?G(-, F) = S Y I,„ for all F e fi, in the sense that 

£ (-G*(-, F)/ + A^DpG jk (; Y)D a cf, i ) = /(F), € C(fi) ffl 

ii) G( , F) e V°''(fi \ g(^^)) for all F € fi and R > 0, and G( , F) vanishes on Sfi. 

iii) For any / = (/',... , / m ) r e C~(fi), the function u given by 

h(X):= f G(F,X)/(F)rfF 
Jn 

belongs to V?' (fi) and satisfies '«5f« = / in the sense of ( 12.51 ). 

Similarly, we say that an m x m matrix valued function G(X, F) = G(f, x, s,y) is a Green's 
function of in fi if it satisfies the following properties: 

i) G(-, F) e W°'l c (Q) and f JSfG(-, F) = 5 r / m for all F € fi, in the sense that 
(g*0, YW, + A^DpGjk{; Y)D a cf> i ) = /(F), € C c °°(fi) m . 



X 



n 

ii) GO, F) e V° a (fi \ g(F,^)) for all F 6 fi and R > 0, and G(-, F) vanishes on Sfi. 

iii) For any / = (/\ ... , / m ) r 6 C™(fi), the function u given by 

k(X):= f G(F,X)/(F)rfF 
Jn 

belongs to V^fi) and satisfies Jz?h = f in the sense of ( |2.4| >. 
We remark that part iii) of the above definitions combined with the uniqueness of weak 



solutions of '«Sfn = / and «5f« = / in V?' J (fi) for any / e C"(fi) gives uniqueness of 



Green's functions; see |4] §3.6] and (3). 



3. Main results 

The following condition (IH) means that weak solutions of 5£u = and '«5f« = enjoy 
interior Holder continuity estimates with an exponent hq. It is not hard to see that this 
condition is equivalent to saying that the operator Jzf and its adjoint LSf satisfy the property 
(PH) in 0; see 12 Lemma 8.2] for the proof. 

Condition (IH). There exist constants fio e (0, 1], R c e (0, oo], and Co > such that for 
all X € fi the following holds: 



GREEN'S FUNCTIONS IN TIME- VARYING DOMAINS 



7 



i) If u is a weak solution of Jz? u = in /?), where R <R C Ad (X), then we have 

[«U/2, W ;e_(X,R/2)<C «^(f |H| 2 ) . 

ii) If h is a weak solution of 'Jz?h = in Q+(X, R), where R < R c A d + (X), then we have 

[«U/2^;G + (X,R/2,<Co J R-"°(-f |M| 2 ) . 

\Jq + (x,r) I 

By assuming the condition (IH), we construct the Green's function of S£ in time-varying 
H\ domains and the domains above time-varying H\ graph in W l+l . 

Theorem 3.1. Let Q. be a time-varying H\ (graph) domain in R' /+1 . Assume the condition 
(IH). Then there exists a unique Green's function G(X,Y) = G(t,x,s,y) of j£f in Q. We 
have G( , Y) e C^ 2 ' m (£l \ {Y))for allY&Q. and 

(3.2) G(-,Y) = on £l(-oo,i). 

Also, there exists a unique Green 's function G(X, Y) of'J£ in Q, which satisfies 

(3.3) G(-,Y) = on Q(i,oo) 

a«(i G(-, F) G Cl^'^iQ. \ {Y})for all Y € Q. In addition, we have the following the identity 

(3.4) G(X, Y) := G{YX) T , MX, Y e Q, X + Y. 
Moreover, for any tff e L2(to(so)) m , the function u(t, x) given by 

(3.5) u(t,x)= ( G(t,x,s Q ,y)iff (y)dy, VX = (t, x) e Q(s , °°), 

Jw(.! ) 

is a unique weak solution of the problem 

(3.6) = in Q.(sq,oo), u — on SQ.(sq,oo), u — if/ on o>(sq) 
and if iff is continuous at xq e o>(sq) in addition, then 

(3.7) Urn | G(t,x, s ,y)tf/ (y)dy = ifr (x ). 



(/,.t)->(i ,-to) 

Xen(sa,co) 



w(s ) 



Furthermore, the following estimates hold for G, where we denote d' Y — d(Y) A R c : 

+ild (Sl\Q(Y,R)) + H G ('' Y )\\V2(CI\Q(Y,R)) ^ NR ^ f° r al1 R < d Y md Y G °' 

ii) HOC-, F)IL p( n KR]) < NR- d+( - d+2 ^ /or a// r < 4, F e fi, andp 6 [1, d -f). 
Hi) \{XeD.: \G(X, Y)\ > t}\ < NT- (d+2)/d for all t > (d' Y /2)- d and Yea 

iv) \\DG(; Y)\\ Lp (n [m ) < NR- d - l+{d+2 »P for all r < d' Y , Y e Q, and p e [1, £f ) . 

v) \{X e Q.: \D X G{X, Y)\ > t}\ < jV T -W+2)/(rf+i) j or all T > ^' Y l2y d and feO. 

vi) \G(X, Y)\ < C\X -Y\~J, whenever < \X - Y\ &> < d' Y /2 and XJefl. 

vii) \G(X,Y)-G(X',Y)\ < C\X-X't^\X-Y\~ d ~^ whenever2\X-X'\a? < \X-Y\& <d' Y /2 
andX,X',Y e il 

In the above, N — N(d, m, v, po, Co) and N depends on p as well in ii) and iv). The estimates 
i) - vii) are also valid for the Green 's function G of the adjoint operator 'Jzf in Q. 

Remark 3.8. In the condition (IH), the constant R c is interchangeable with aR c for any 
fixed a € (0, 00), possibly at the cost of increasing the constant Co- Also, the condition (IH) 
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implies that if u is a weak solution of _£? u — in Q-(Xo, R) with R < d (F) A R c , then we 
have the Loo estimate 

(3-9) H»llL„(e-(W))<ivff |h| 2 ) , 

\JQ-(Xo,R) I 

where N = N(d, m, v,po, Co) > 0. Moreover, u satisfy 

l|H|k( G _ ( x„,-)) < N(R - r)- (d+2) IP\\u\\ Lp{Q _ {Xa , R)) , Vr < R, Vp> 0, 
where N = N(d, m, v, po, Co, p) > 0. See J4] Lemma 2.6] for the proof. 



Remark 3.10. In Theorem 13. II we also have the following estimates, which follow from 
the identity ( 13.4-b and the estimates i) - vi) for G(-,X): 

V \M X > Oil W n\e ( x, fi)) + \\G(x, -^(nmrn) * NR ~ dl2 for a11 R < d * md x e a 

ii) \\G(X, ■)\\ LpmxM < NR- d+{d+2) IP for all r < d' x , X e O, and/? e [1, ^). 

iiij |{F e £2: |G(X, F)| > t}| < Nr- (d+2)/d for all t > (d' x /2T d andX € £1 

/v) ||DG(Z, -)lli,(PPWD ^ NR- d - 1+( - d+x »P for all r<4Jefl, and p e [1, ffi) . 

vj |{F € £2: |D V G(X, F)| > t}\ < N T -V+md+i) for all T > (d' x l2T d and lefl. 

vj'J |G(X, F)| < C\X - Y\~ d whenever < \X - Y\&> < d' x /2 and X, Y e £1 

vz7j |G(X,F)-G(X,F')| < C|F-F'|*|X-F|^~ W whenever 2\Y-Y'\g, < \X-Y\& < d' x /2 
and X, F, F' e £1 

In particular, |G(X, F)| < - F|^ whenever < |X - Y\& < \(d(X) V d(F)) A R c . 

The following condition (LB) is used to obtain a global Gaussian bound for the Green's 
function G(X, Y) in a time-varying H\ domain Q. c R d+1 . 

Condition (LB). There exist constants R max e (0, oo] and A^o > so that for all X e O and 
< R < R max , the following holds. 

i) If u is a weak solution of Jzfu = in Cl [X, R] vanishing on 3^Q.-[X, R], then we have 

ll«lk.<n_PW?D < ^" (2+rf)/2 ll«llL 2( n_[x,R])- 

ii) If u is a weak solution of 'S£ u — in Q+ [X, /?] vanishing on ^£2+ [X, 7?] , then we have 

IImIIl^rk/2]) < N R- {2+dy2 \\u\\ L2(n+[XR]) . 

Theorem 3.11. Let Qbe a time-varying H\ (graph) domain in R . Assume the condition 
(LB) as well as the condition (IH). Then the Green 's function G(X, Y) of ' ££ in Q exists and 
satisfies the conclusions ofTheorem \3.1\ Moreover, for all X — (t,x) and Y — (s,y) in £2 
with t > s, we have 

(3.12) \G(t, x, s,y)\ <N{(t- s) A RlJ"" 2 exp [-k\x - y\ 2 /(t - s)} , 

where N — N(d, m, v, A^o) and k — k(v) > 0. 

Remark 3.13. In the condition (LB), the constant R max is interchangeable with aR max for 
any fixed a e (0, oo), possibly at the cost of increasing the constant No. In Theorem l3.1 11 
the estimate ( 13.121 ) implies, via straightforward computation, that 

(3.14) \G(X,Y)\<N\X-Y\~ d , if < |f - *|< R 2 ^, 
where N = N(d, m, v, Nq). Then, similar to Lemma l5T8l below. one can show 

(3.15) ||G(, Y)\\ LiMmm) + ||G(, Y)\\ V2inmM < NR~ dl2 , VR e (0,R max ), 
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where N = N(d,m,v,No). Moreover, using ( 13.151 and proceeding as in Section 4.2], 
one can show that G satisfies the estimates ii) - vi) in Theorem 13. II with d' Y replaced by 
Rmax- Also, it is clear that the estimate ( 11. 2t in the introduction follows from Theorem l3.1 II 



In order to derive the estimate dl.3l ) in the introduction, we introduce the following 
condition (LH) which, loosely speaking, says that weak solutions of ££ u = and u = 
vanishing on E c 5Q are locally Holder continuous up to E with exponent //o- 

Condition (LH). There exist fi Q e (0,1], R max e (0, oo], and N\ > so that for all X e Q 
and < R < R max , the following holds. 

i) If u is a weak solution of — in Cl-[X, R] vanishing on £? > Cl-\X,R~\, then we have 

[H]^/2,^;qjx,r/2) < NiR^° i -f \uf\ , where u =Xn_[xji]U. 

\Jq-(x,R) I 

ii) If u is a weak solution of 'S£ u — in Q + [X, R] vanishing on [X, R] , then we have 

[«W2 rt ;e t (W2)<WiF''»f-f |h| 2 ) , where u =*n + [XR]«- 

\JQ + (X,R) I 

It is easy to see that the condition (LH) implies the condition (LB); see Lemma loTTl in 
Appendix for the proof. Also, it is obvious that the condition (LH) implies the condition 
(IH). Therefore if the condition (LH) is satisfied, then there exists the Green's function of 
S£ and it satisfies the conclusions of Theorems 13 . 1 1 and 13 . 1 1 1 The following theorem says 
that in fact, in such a case, a better estimate for the Green's function is available. 

Theorem 3.16. Let £lbe a time-varying Hi (graph) domain in R d+1 , Assume the condition 
(LH). Then the Green's function G(X, Y) of J*f in Q. exists and satisfies the conclusions of 
Theorem \3.1\ Moreover, for all X — (t, x) and Y — (y, s) in Q. with t > s, we have 

(3.17) \G(t, x, s,y)\ < N6(X, Yf%t - s) A R 2 max r dl1 exp {-k\x - y\ 2 /(t - s)} , 
where N — N(d, m, v, /j.q, N\ ) and k — k(v) > and we used the notation 

( 3.18) S(X, F) = (l a ^ , Via d+m 



Rmax A IX - Y\» ) \ R ma x A |X - Y\ & 

Remark 3.19. In the condition (LH), the constant Rmax is interchangeable with aR max for 
any fixed a e (0, oo), possibly at the cost of increasing the constant Ni . Also, we note that 
the estimate ( 11.3b in the introduction follows from Theorem l3.16l if Q be a time-varying Hi 
domain or a time-varying Hi graph domain with R max = oo. 

Remark 3.20. In Theorem l3.16l we also have the estimate 

NS(X,Yy> (IX-X'l&Y ( *|x-y| 2 l 

G(X,Y)-G(X',Y) < — —{■ — exp<^ — 1 —\ 

{{t-s)KRl ax } d!2 \\X-Y\&) V \ t-s J 

whenever 2|X -X'\&> < |X- Y\& and t > s. It follows from ( 13.171 ) and the condition (LH). 



4. Some Applications of Main Results 

4.1. Scalar case. In the scalar case (i.e., m = 1), both conditions (LB) and (IH) are sat- 
isfied with R c = Rmax = oo and A^o = No(d,v); see e.g., |22, Chapter VI]. Also, in the 
scalar case, the Green's function is a nonnegative scalar function. Therefore, the following 
corollary is an immediate consequence of Theorem l3.1 II 
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Corollary 4.1. Let CI be a time-varying H\ (graph) domain in M. d+l . If m — 1, then the 
Green 's function G(X, Y) of Jz? in CI exists and satisfies the conclusions of Theorem 13.71 
with d' Y replaced by R a . Moreover, for all X — (t, x) and Y — (y, s) in CI with t > s, we have 



G(t,x, s,y) < N(t- sy a,A exp{-K\x-y\y(t- s)} , 

where N — N(d, v) and k — k(v) are universal constants independent of CI. 

In fact, in the scalar case, a better estimate is available near the boundary. Let CI be 
a time-varying Hi (graph) domain in R d . By using the results in l22l §VI.8], one can 
show that in the case when m - 1, the condition (LH) is satisfied in Q. Moreover, in the 
case when Q is a time-varying Hi graph domain, then the condition (LH) is satisfied with 
Rmax = 00 ■ Also, in that case, there exists TV = N(M) > 1 such that 

(4.2) 1 < d-(X)/d(X), d + (X)/d(X) < N, Vie Q. 

Therefore, the following corollaries are immediate consequences of Theorem l3.16l 

Corollary 4.3. Assume that m — 1 and let G(X, Y) be the Green 's function of Jz? in CI, 
where CI is a time-varying Hi domain in M. d+l . Let 5(X, Y) be as defined in ( 13.18b with 
Rmax — Ra- Then, for all X — (t, x) and Y — (y, s) in CI with t > s, we have 

G(t, x, s,y) < N5{X, Yf"{(t - s) A R 2 a }- d/2 exp {-k\x - y\ 2 /(t - s)) , 

where N — N(d, v) and k — k(v) are positive constants independent of CI. 

Corollary 4.4. Assume that m — 1 and let G(X, Y) be the Green 's function of I£ in Q, 
where Q is a time-varying Hi graph domain in R d+ 1 . Then, for all X — (t,x) and Y — (y, s) 
in Q with t > s, we have 

I d(X) Y°( d(Y) Y° 1 f K\x-y\ 2 

G(t,x,s,y)<N 1A -— 1A -— TTrexp^ — ! — 

y ' \ \X-Y\&) \ \X-Y\&) (t-s)*' 2 P \ t-s 

where N — N(d, v, M) and k — k(v) are positive constants. 

4.2. Leo -perturbation of diagonal systems. Let a al3 (X) be scalar functions satisfying 

d 

(4.5) a^(X)^ a > vo|^| 2 , V^elT'; £ \a^(X)\ 2 < vf, 

for all X € R d+l with some constant vq e (0, 1]. Let Q be a time-varying Hi (graph) domain 
in R l/+1 . Let A a ? be the coefficients of the operator J£. We denote 

( „ , ^ 1/2 



(4.6) #= sup \ YMjW-^w^ii 



\i,j=laj}=\ 



where 6ij is the Kronecker delta symbol. By Lemma 15781 there exists Sq = So(d, vq,M) 
such that if S < So, then the condition (LH) is satisfied with po = po(d, vq, M), R max = R a , 
and N\ = N\{d,m,V(j,M). Therefore, the following corollaries are another easy conse- 
quences of Theorem |3.16| 

Corollary 4.7. Let CI be a time-varying Hi domain in R d+1 and let S(X, Y) be as in ( 13.18| > 
with Rmax — Ra- There exists So — So(d, vo,M) such that if S < So, then the Green's 
function G(X, Y) of J?f in CI exists and satisfies the conclusions of Theorem 13.71 with d' Y 
replaced by R a . Moreover, for all X — (t, x) and Y — (y, s) in CI with t > s, we have 

\G(t, x, s,y)\ < NS(X, Yf°{(t - s) A R 2 a Y dl2 exp{- K \x - y\ 2 /(t - s)}, 

where N, po, and k are constants depending on d, m, vq, and M. 
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Corollary 4.8. Let CI be a time-varying Hi graph domain in R . There exists £§ — 
$o(d,vo, M) such that if S < Sq, then the Green's function G(X, Y) of _Sf in Q exists and 
satisfies the conclusions ofTheorem \3.1\ with d' Y replaced by oo. Moreover, for all X — (t, x) 
and Y — (y, s) in Q, with t > s, we have 

I d(X) Y" 1 / d(Y) Y"° 1 f K\x-y\ 2 ) 
\G(t,x,s,y)\<N 1 A — — 1A -— Tj^expl — 1 —}, 

where N,po, and k are constants depending on d, m, vo, and M. 

4.3. Systems with VMO A coefficients. For a measurable function / = f(X) = f(t,x) 
defined on R d+l , we set for p > 

(o p (f) := sup sup + j- \f(y,s)-f v (s)\dyds; f x , r (s) = -f f(s,-). 

XeR' i+1 r<p Jt-r 2 JB(x,r) J B(x,r) 

We say that / belongs to VMO. t if lim p ^ <Q p (f) = 0. Note that VMO A - is a strictly larger 
class than the classical VMO space. In particular, VMO A contains all functions uniformly 
continuous in x and measurable in t ; see Ifl9l . 

By J4] Lemma 2.3], we find that if the coefficients of Jz? belong to VMO A , then the 
condition (IH) is satisfied with parameters ^o. No, and R c depending on m p (A a ^) as well as 
on d, to, v. Therefore, we have the following corollary of Theorem |3.1| 

Corollary 4.9. Let Q be a time-varying H\ (graph) domain in R d+1 . If the coefficients of ££ 
belong to VMO A -, then the Green' s function of J*f exists in Q. and satisfies the conclusions 
ofTheorem \3. 1 1 with some R c > 0. 

Remark 4.10. In Corollary 14.91 instead of assuming A a ^ € VMO*, one may assume that 
iQ p {A a P) is sufficiently small for some p > 0. Also, if Q is a time-varying domain satisfying 
the hypothesis of Section 2.3 with / = f(X') = f(t,x') 6 Hi +a (R d ) for some a > 0, then 
one can show that the condition (LH) is satisfied with the parameters po, N\, and R max < 00 
depending on d, to, v, and <s> p (A a ^); see ll22l for the definition of the space H\+ a . Therefore, 
in that case, for all X = (t, x) and Y — (y, s) in £1 with t > s, we have 

\G(t, x, s,y)\ < NS(X, YT\{t - s) A R^r" 2 exp{-*|x - y\ 2 /(t - s)}, 

where 6(X, Y) is as in (I3T81 . 

5. Proofs of Main Theorems 

5.1. Proof of Theorem 13.11 By following |4), we shall first construct the "averaged" 
Green's function of _S? in £1 Notice that we have dD. = SPCl = <?£l = SQ.. The fol- 
lowing lemma is used for the construction of the averaged Green's function, which follows 
essentially from Brown et al. [3] and an embedding theorem in ||20l §11.3], which is also 
valid for functions in P2' (Q(£q, fi))- We remark that the function space Vy 1 ^) coincides 
with the function space Vo(^) used in J3). 

Lemma 5.1. For g e Q"(Q(fo, h)) m and if/ Q € L2(co(to)), there exists a unique weak 
solution v e V°''(Q(fo, h)) of the problem 

v = g in Q(fo, fi), v = on SCl(to,ti), v = i/f on w(fo)- 

Moreover, we have the following energy inequality for the weak solution v: 

(5.2) l|v||y 2 (o ((b ,„)) <N^\g\\ k2dUMMm , 0ttl)) + ||^olli2(w&o))); N = N(d,m,v). 
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Similarly, for f € C£°(Q(fo, t\)) m and \\i x e L2(w(fi)), there exists a unique weak solution 
u e V°' l (Q.(t , ti) of the problem 

'Jz?h = / in Q(fo,fi), u = on 5T2(fo,fi), u -ij/ x on a>(ti) 

and u satisfies the following energy inequality: 

IMIvynte,*!)) < N (\\f\\ L{2d+mM)mo , tl )) + llWktetfti))); # = N(d,m,v). 
Furthermore, we have the identity 

(5.3) I f-vdX+ I v \fi x dx = I u gdX+ I utff dx. 

Let us fix a function € C~(R rf+1 ) such that <E> is supported in Q-(0, 1), < < 2, and 
0=1. Let Y = (s,y) e Q. be fixed but arbitrary. For < s < d{Y), we define 

£ (X) = £ (f, x) = £- d " 2 0((t - s)le\ (x - y)/s). 

Fix to € (-oo, 5 - s 2 ) and let v = v E jk be a unique weak solution of the problem 

Jzfv — E ek in fi(fo,°°), v — on SiXfo, 00 ), v = on <w(?o), 

where e k is the £-th unit vector in R m . By the uniqueness, we find that v does not depend 
on the particular choice of to and we may extend v to the entire Q by setting 

(5.4) v = v eXk =0 on Q(-oo, s - e 2 ). 
Then v e V®' 1 ^) an d satisfies for all t > s the identity 

\ v'ftdx- \ v i 4dX+ \ A^Dpv 1 D a ^ i dX= \ d> s <p k dX 

for all e C™(Cl) m . We define the averaged Green's function G e (-, Y) = (G B jk (; Y))"J k=l of 
the operator Jz? in Q by 

G ^ ( ,y) = v i = v ^. 
Notice that by Lemma lBTTI and an embedding theorem (see 11201 §11.3]) we obtain 

(5.5) \\G E (; Y)\\ L2+4lAn) < N\\G%, Y)\\vm < N\\®s\\ k ^ )KM) (.n) < N £- d( ' U2) l (2d+4) . 

Next, for any given / € C^°(Q)" ! , fix t\ such that / = on Q(fi, oo) and let u be a unique 
weak solution of the problem 

'•ifii = / in Q(-oo,f[), ii — on SQ.{-co,t\), u — on u{t\), 

Again, by the uniqueness we may extend u to the entire Q. by setting u = on 0(fj, oo). 
Then, u e V%' l (Cl) and satisfies for all r the identity 

| u i 4> i dx+ \ u i <p t dX+ \ A"fD p v i D a <t> i dX= \ f'fidX 

for all 6 C£°(Q). Also, similar to (15.5b . we have 

1 1" 1 1^2+4/,/ (Q) < Nil/Ill 

Now, let Xo € O and 7? < d(Xo) A be fixed but arbitrary, and assume that / is supported 
in Q + (Xq,R) c SI. By using the condition (IH) and following the same argument as in |4] 
Section 3.2], we obtain 

(5.6) I|h|L, (G+ (Xo,*/4)) < NR 2 -^ 2)lp \\f\\L„(Q +( x ,R)), Vp>{d + 2)12. 
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If Q-(Y, s) c e+(X , R/4), then (O together with $5M yields 

f G £ (-,T)/ < f O e |h| < ^ 2 - w+2)/ ni/ll VG+ (x„,«)), -ip>{d + 2)l2. 

JQ + (X ,R) JQ-(Y,e) 

By duality, it follows that if e) c <2+(Z , J?/4), then we have 

\\G B (; Y)\K( Q+ (.Ko,R)) < NR- d+(d+ W, Vq e [1, (d + 2)1 d). 

Then, by following the proof of [4, Lemma 3.2], we conclude 

(5.7) \G B (X, Y)\ < N\X - Y\J, < \\X - Y\& if \X - Y\& < \(d(Y) A R c ). 

The following lemma is a consequence of the energy inequality of Brown et al. |3), the 
above estimate ( 15.7b . and an embedding theorem in 11201 §11.3]. 

Lemma 5.8. For R < \(d(Y) A R c ), let £ € C™(Q(Y,R)) be a cut-off function such that 
Q<(<\and£=lon Q(Y,R/2). Then, for all e > we have 

lid - 0G e {-, Y)\Ww) ^ MIIATIlL + Il6ll0 1/2 « 1 " d/2 - 

In particular, for all s > and R < j(d(Y) A R c ), we have 

\\G%,Y)\\ V2(nmzR)) <NR- d/2 . 

The following lemma is an analogue of |4] Lemma 6.1] in time-varying H l domains, 
the proof of which is essentially the same. 

Lemma 5.9. Let {uk}^ be a sequence in V2(Q). 7/'sup i: ||Mj : ||y 2 (£x ) < N < oo, then there 
exists a subsequence {u kj }J =l £ {uk}'^L l and u 6 V2(Q) with ||w||y 2 (n) < N such that u/ ii — k u 
weakly in W^' l {Q.{k), t\))for all — oo < f < t\ < oo. Moreover, if all u^ vanish on SO., then 
u also vanishes on S £1 

The above two lemmas contain all ingredients needed for the construction of a Green's 
function. By following the argument in O Section 4.2] verbatim, we construct the Green's 
function G( , Y) from G e (-, Y), and it is readily seen that G(-, Y) e Cf o ° c /2,Wl (Q \ {Y}) satisfies 
( 13. 2t as well as the estimates i) - vi). The estimate vii) does not appear explicitly in |4| but 
it easily follows from the estimates vi) and the condition (IH); see |[T4l §3.6]. 

Also, fix a function *F e C™(R d+1 ) such that ¥ is supported in g+(0, 1), < *F < 2, and 
jL, +1 *F = 1. For < e < d(Y), where Y — (s,y) e D. be fixed but arbitrary, we set 

V E (X) = ^s(t, x) = s- d - 2 V((t - s)ls 2 , (x - y)ls). 

Fix 1 1 € (t + s 2 , oo) and let w = w e xk be a unique weak solution of the problem 

'Jz5f = W s e k in Q(-oo,r,), w = on Sfit-oo,^), w = on uXf\), 

Then, as before, we may extend w to the entire Q. by setting w = on 0(f + e 2 , oo) so that 
w belongs to Vi'^Q) and satisfies for all r < t the identity 

(5.10) \ w i (j) i dx+ \ w i (j) i ,dX+ \ AfDpw j D a (l) i dX= \ V e <f> k dX 
Jw(t) JntT.M) Jn(T,oo) /; Je + (x,s) 

for all 4> e C~(Q)"'. We define the averaged Green's function G £ (-, Y) = (G s jk (; Y))J k=1 of 
the adjoint operator '££ in Q by 
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Then by a similar argument, we construct a Green's function G(-, Y) from G e (-, Y), which 
belongs to Cf° ^(Q \ \Y}) and satisfies (I3.3l i and the estimates /J - vi). Moreover, by 
following [4 Lemma 3.5], we obtain the identity (13.4I >. 

Next, we shall prove the identity (13.5b . Let ^ e L2(w(so))'" be given and let u be a 
unique weak solution of the problem ( 13.6b . Fix X = (f, x) e C2( so, °°) and let w = w Et x,k be 
as constructed above. By Lemma lBTTl for s sufficiently small, we have 

(5.11) f y B u k dY= f w E , x .k-%dy. 

jQ + (X,e) JoKso) 

For if/ € C"(w(so)), it can be easily seen that (see |4] Section 3.5]) 
lim I w B wi/r dy= G(-,X)e k ■ t/f dy, 



Since the condition (IH) implies that u is continuous at X, by taking the limit e — > in 
( 15. lit and using ( 13.4b , we obtain 

u*(X) = I G,i.(so,;y,f,-*)(/'oO ; M)' = I G k i(t,x,s ,y)^ (y)dy. 

Jii)(s ) Jio(s ) 

We have thus derived (13.5b under an assumption that ^ € C^(w(so))- For t// Q e L2(w(so)) m , 
let {*^j}J =l be a sequence in C^(a>(so)) m such that ^ — > ^ in L2(w(so))- Let My be a 
unique weak solution of the problem (13.6b with ^ replaced by Then by Lemma lBTTl 
we find that lim^ooHitj - u\\v 2 (n(s ,t)) = and by the condition (IH) and ( 13.9b we have 
limj^cc\Uj(X) - u(X)\ = 0. On the other hand, by the estimate i) applied to G(-,X) together 
with the identity ( 13.4b , we find that \\G(t, x, sq, •)IIl 1 (w(.! )) < °°> ar, d thus we get 

lim I G(t,x,s ,y)fj(y)dy= I G(t,x, s ,y)tff (y)dy. 

This completes the proof of (I3.51 l. Similarly, for ijr l e L,2{u{t\)) m , let h be a unique weak 
solution of the problem 

'J? v = in Q(-oo,^), v = on SQ(-oo,^), v = >fr l on w(/i). 

Then as above, v has the following representation: 

v(s,y) — I G(s,y,ti,x)>ff l (x)dx. 

JaXti) 

It only remains us to prove ( 13.7b . We proceed similar to J4] Section 4.4]. The following 
lemma is another simple consequence of Brown et al. Q. 



Lemma 5.12. Let rj — rj(x) e C (R ) be a bounded nonnegative function. Assume that 

■o,: 

2 



u G V^' 1 (Q.(sq, oo)) is the weak solution of the problem (13.6b and define 



/(f) = - f /7(x)|«(f, x)\ 2 dx, te(s ,oo). 
77ie« 7(f) is absolutely continuous and satisfies a.e. t > sq the identity 



I'(t) = - f A^D p u j D a {7j^). 
The following lemmas are key ingredients to prove ( 13.71 ) and adapted from |4|| . 
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Lemma 5.13. Assume that if/ e L2(to(so)) m is supported in a closed set F c ll>(sq) and let 
u be the weak solution of the problem ( 13.61 ). Then, we have 

(5.14) f\u{t,x)\ 2 dx <e- ydisl{E ' F)2/{ '- So) f\ij/ Q {x)\ 2 dx, V£cw((), 

Je Jf 

where dist(£, F) = inf{|x — y\: x e E, y e F] and y = y(v) > 0. 

Proof. We may assume that dist(£, F) > 0; otherwise ( 15.141 ) is an immediate consequence 
of the energy inequality (15.2b . Let (p = <p(x) be a bounded C 1 function on W 1 satisfying 
\D<f>\ < K for some K > to be fixed later. Define 

1(f) = f e 20W |H(f,-*)| 2 ^, ?>io- 

By Lemma [5 .121 we find that /'(f) satisfies for a.e. f > so 

I'(t) = -2 f e^A^Dp^D^dx-A f e^A^Dpu'D^u 1 dx 
J(jj(t) Jiit(t) 

<-2v f e 2<f, \Du\ 2 dx + 4(K/v) f e^\Du\e^\u\ dx 

J(o(t) J0}(t) 



off) J (off) 



(2/v i )K 2 f e 2( *|«| 2 = (2/v")K%t). 

Ja>(t) 



The above differential inequality yields 

(5.15) 7(f) < e (2/v3) ^ ( '- So, lkVollL(F)' Vf ^<>. 

Notice that by a standard approximation, we may assume that is a bounded Lipschitz 
continuous function satisfying \D<f>\ < K a.e. Since F is a closed set, the function 

dist(x, F) = inf{|x - y| : yeF) 

is a Lipschitz function on R d with Lipschitz constant 1 and dist(£, F) = inf x€ e dist(x, F). 
Therefore, if we set <p(x) - K(dist(x, F) A dist(£, F)), then by (15.151 ), we get 

j \u(t, x)\ 2 dx < exp{(2/v 3 )K 2 (t - s ) - 2K dist(E,F)] j \ifi (x)\ 2 dx. 

The lemma follows if we set K = dist(£, F)/{(2/v 3 )(t - s )}. ■ 

Lemma 5.16. Let u be the weak solution of the problem ( 13.6b , where if/ e L^oXso)) and 
has a compact support in lj(sq). Denote 

(5.17) g = g(x) = dist(jc, dcu(so)) A R c , x e u>(so) 

Then for all x e oj(soX we have 

\u(t,x)\ < A^||^ || Lk , Mso) ) whenever < t- s < (1 A M~ 2 )g 2 (x)/4, 
where N — N(d, m, v, /j.q, Co) > 0. 

Proof. For x e a>(so), set r = g(x)/2 and 6 = (r/M) 2 so that 

{so -S,so + S)x B(x, r) cc O. 
For any t satisfying < t - so < 5 A r 2 = (1 A M~ 2 )g 2 (x)/4, set R = -\Jt - so and denote 
A = B(x,R); A k = {y e o)(s ): 2 k ~ l R < \y - x\ < 2 k R), A: = 1,2, 
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Since i// is compactly supported in cj(so), we have i// = J^ =0 XA k ^o f° r some Icq < oo. For 
k = 0, 1, . . . , ko, we define 



u k (t, x) = I 

JA k 



G(t, x, s ,y)if> (y)dy. 



Then, it follows from (13.5b that u = «a and that each u k is the weak solution of the 
problem ( 13.6b with^,^ in place of ^ . We apply Lemma lBTTI to u k with £ = B(x,R) and 
F = for £ = 1 , 2, . . ., to obtain that 

f |«*(s,y)| 2 rfy < ^- r(2 " , - I, 2 fa V||^ ||^ (w(!o)) , V, e (s , 0- 

JB(x,R) 



ll^olk < A^ll^olk 



Therefore, by the condition (IH) and ( 13.91 ). we get 

k I oo ^ 

< £|«,(f,x)| < AT 1 + ^ e -r(2*--i) 2 W2 

<r=0 V k=\ t 

The lemma is proved. 

Lemma 5.18. Let r\ e C™(B(xo,4r)) be a function satisfying 

(5.19) < i] < 1, ^il /« B(x ,2r), and \Dr]\ < 4/r, 

where xq € a>(so) a«of r < p(xo)/5, where g is as defined in ( 15.17b . 77ie«, we /zave 



lim r 



G(t, x, s , y)r](y) dy = I,„, Vx e B(x , r). 

so) 



Proof. By taking £ — > in ( 15.10I ) and arguing as in the proof of [4', Lemma 4.3], we find 
that the following identity holds for all r < t: 

(5.20) <l> k (X) = <p k {t,x)= f G ik (T,y,t,x)4> i (T,y)dy+ f G ik (Y, X)^(F) dY 

+ f A; ! fD y/i G 7 ,(T,XW'(T) £ /F, e C f °°(C2) m , 

Jfi(T.f) 

where we have used ( 13. 31 ). Let £ = £(s) be a smooth function on R such that 

< £ < 1, £(s) = 1 for |s - \ I < <5, and £(s) = for \s - s \ > 25, 

where 5 is chosen so small that 

(so - 25, s + 25) x B(x , 4r) cc Q. 

Notice that we may take 0(T) = 0(s,y) = £(s)ri(y)et in ( 15.20b . Setting t = so and assuming 
that \t - s | < 5 in ( 15.20b . we obtain by ( f3~4T > that 

(5.21) S kl = f G k ,(t,x,s ,y)T](y)dy+ f A^D, fi G Jk {Y,X)D a Ti(y)dY =: / + //. 

Then for all X = (f, x) such that x e Z?(xo, r) and |? - sol < 5 A r 2 , we estimate 77 as follows 
by using the hypothesis ( |5.19b . Holder's inequality, and the estimate i) for G(-,X): 

\II\ < Nr djl -\t - s ) 1/2 ( f \D v G(Y,X)\ 2 dY] < Cr-\t - s ) 1/2 . 

\jSl(s ,t)\Q(X,r) I 

Therefore, the lemma follows by taking the limit t to so in ( 15.211 ). ■ 
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We are ready to prove d3.71 i. Let e L2(u>{so)) m and assume that if/ is continuous at 
xo e u>(so). Let u be the weak solution of the problem ( 13 .6b - For any e > given, choose 
r < q(xq)/5, where g is as defined in (15.17b . such that 

\^q(x) - ^o( x o)l < s/2N for all x satisfying \x - xq\ < Ar, 

where N is the constant that appears in Lemma 15. 161 Let r\ be given as in Lemma 13. 18 l and 
let «o, u e , and h^, respectively, be the weak solution of the problem ( 13. 6b with T]i// (xo), 
T](if/ - >p (xo)), and (1 - rf)if/ in place of if/ . By the uniqueness, we have u = uq + u s + u„, 
and by the formula (13.5b . «o is represented by 



Uo(t, x) = I 



(5.22) ua(t,x) = \\ G{t,x,so,y)Ti(y)dy\i// {x Q ). 

^w(so) / 

Let 6 > be chosen so that 

(s - 5, s + 6) x B(x , Ar) cc Q. 



For any s satisfying < s - sq < S, we set £ = B(xo, r) c w(i) and F = u>(so) \ B(xq, 2r) 
in Lemma l5.13l to get 



JB(.v ,r) 



Therefore, for all f satisfying < f - sq < 6 A r 2 , we set R = y? - so/4 in ( 13.9b to get 

(5.23) |n»(t,x)| < Mr^V* r/ * )2 ||^ ||^ M , Vx£B(x ,R). 
Finally, we estimate u e by using Lemma IB. 161 

(5.24) \u E (t,x)\ < s/2 whenever < /- s < (1 A M~ 2 )g 2 (x)/A. 

Combining (15.22b , ( 15.23b . and (15.24b . we see that if t — sq is chosen sufficiently small, 
then there exists § > such that for all x e B(xq, §) we have \u(t, x) - ifro(xo)\ < £■ This 
completes the proof. ■ 

5.2. Proof of Theoreni l3.11l By Theorem l3.ll the condition (IH) implies existence of the 
Green's function G(X, Y) of Jz? in Q. Therefore, all the conclusions of Theorem 13. II are 
satisfied. Let <p be a bounded Lipschitz function on R d satisfying \D<p\ < K a.e. for some 
K > to be chosen later. For any / € La((o(s)) m , let u be a unique weak solution of the 
problem 

(5.25) = in Q(s,oo), « = on 50(5,00), u = e^f on a>(s). 
For t > s, we define the operator : L2((o(sJ) m — > L2(a>(t)) m by 

Notice that by the representation formula (13.5b . we have 

(5.26) = f Git,x,s,y)e-^f{y)dy. 
For f > i, we define 

/(o = ii^«(^oii^ (0) = ii^,/iiL ( . W )- 

Then, as in the proof of Lemma 15.131 we find that / is absolutely continuous and satisfies 
for a.e. t > s, the differential inequality 

/'(f) < (2/v 3 )K 2 I(t). 
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The above inequality with the initial condition I(s) = \\f\\'[^ oj(s)) yields 

7(0 < ^^ 2 ^>||/Hi 2Mj)) , W > s. 
We have thus shown that for all / 6 L2(u){s)) m , the operator P^ t satisfies 

(5.27) Wi+JKWS) e^^WfWuaw Vf > s < 

where § = v~ 3 . We set R = yjt - s A R max and use the condition (LB) to estimate 
e-^ x) \PU t f(x)\ 2 = \u(t,x)\ 2 = \u(X)\ 2 

<N 2 R <d+2) \ \u(Y)\ 2 dY 

= N 2 R -id + T> C' r e - 2 ^\PU T f(y)\ 2 dydr, 

Jl-R 1 Joj(T)riB(x,R) 

Therefore, by using the estimate d5.27t , we get 

\PUf(x)\ 2 <N 2 R- d - 2 f f e 2 ^- 2 ^\PU T f(y)\ 2 dydr 

Jt-R 2 Jaj(T)nB(x,R) 

<N 2 R-"- 2 f f e 2KR \PU T f(y)\ 2 dydT 

<N 2 R- d - 2 e 2KR f e^^Wf^dT 

<N R e II/ILm,))- 
We have thus obtained the following L2(co(s)) — > Loo(w(0) estimate for P^ t : 

(5.28) WPUfWu^W) < N R- d/2 e KR+ ° K2 ^\\f\Uo*s)). 

We also define the operator Qf_, s '■ L2(co(t)) m — > L2(to(s)) m for s < t by setting 

where v is a unique weak solution of the problem 

(5.29) 'Jzfv = in Q(-oo,f), v = on 5Q(-oo, f), v = e*g on Uf). 
By a similar calculation that leads to ( 15.281 ), we obtain 

(5.30) WQUsSUl-M*)) ^ N R- d/2 ^ fi+ ^ 2( '- s) ||g|| L2M ,))- 
It follows from ( 15.251 ), ( 15.291 ), and d531 > in LemmalBTTlthat 

f (Pt^,f)-gdx= f f-{QUs8)dx, V/6L 2 Mi))'", Mg e L 2 (co(t)) m . 

Jw(r) Jbj(s) 

In particular, the above identity holds for all / € C™(aj(s)) m and g e C™(co(t)) m . Therefore, 
by the estimate ( 15.30b and duality, we get 

(5.31) WPUfWma*,)) < N R- d/2 e KR+ » K2{, - s) \\f\\L M V/ e C c °°(<^)) m . 
Now, set r = (s + t)/2 and observe that by uniqueness, we have 

PUf = PU,(PUrf), V/ e C™{gj{s)T. 

Then, by noting that t — r = r — s = (t — s)/2 and R/ V2 < -\Jt - r A R mdx < R, we obtain 
from (EH) and (ED) that 

WpLJWl^ < NR- d e 2 ^ + ^ 2( '- s) ||/|| LlMj) ), V/ e C~{io{s)T; N = 2 d/2 N 2 . 



<h(r) 
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For all x e a>(t) and y e u>{s), the above estimate and (15.261 yield, by duality, that 

(5.32) e^-^Git, x, s, y)\ < NR~ d e 2 ™^^. 

Let (p(z) = K(po(\z - y|), where tpo is defined on [0, oo) by 

\r if r < \x — y\ 

\x-y\ if r>\x-y\. 

Then, <p is a bounded Lipschitz function on R d satisfying \D(f>\ < K a.e. We set 

K= \x -y\l2§{t- s) and £ := \x - y\/ tJF^~s. 

By ( 15.321 ) and the obvious inequality R/ y/t - s < 1, we have 

\G(t, x, s,y)\ < NRT d exp{£/i? - f/4§}. 

Let N = N(&) = N(v) be chosen so that 

exp(£/# - < N exp(-£ 2 /8#), V£ € [0, oo). 

If we set k — 1 /8# = v 3 /8, then we obtain 

\G(t, x, s,y)\ < NR- d exp \-k\x - y\ 2 /(t - s)] 



where N = N(d, m, v, No) > and recall that we set R = y/t - s A R mnx . The proof is 
complete. ■ 

5.3. Proof of Theorem 13.161 Notice that by Lemma [67TI and Theorem |3.1 II for all X = 

(t, x) and Y - (y, s) in O with t > s we have 

(5.33) \G(t, x, s,y)\ < Ci {{t - s) A S^Jf*' 2 exp {-k\x - y\ 2 /(t - s)} , 

where C\ = Ci{n,m,v,^,o,N\). We denote 

*C*,n = (lA„ ™ m ) and M F) = (lA ^ 



^ma.v A \X - Y\gn ) \ R m ax A \X - Y\gn 

so that 6(X, Y) = 6i(X, Y) 6 2 (X, Y). To prove the estimate ( 13.17b . we first claim that 

(5.34) \G(t, x, s,y)\ < NS^X, Yf° {(r — s) A R 2 ^}"" 2 exp [-k\x - y\ 2 /4(t - s)\ , 

where N = N(n, m, v, yUo, N\). The following lemma is a key to prove the above claim. 

Lemma 5.35. For R e (0,R max ) andX e Q such thatd~(X) < R/2, letu be a weak solution 
ofJtfu — in £l-[X, R] vanishing on &Q.-[X,R]. Then, we have 

(5.36) \u(X)\ < Nd-(Xf'>R- d/2 - 1 -'">\\u\\ L2(n _ lx , R]h 
where N = N(d, m, v, fi Q , Ni , M). 

Proof. By the very definition the condition (LH), we have 

(5.37) \u(Y) - u(X)\ < N\Y - X\% B^^^M^^, VF e Q-(X,R/2). 

For any r satisfying d~(X) < r < R/2, there is F e Q_(X,R/2) \ Q such that \X - Y\g? = r. 
By ( 15.371 ) we obtain 

\u{X)\ = \u(X) - u{Y)\ < NrWR-<VZ-i-»°\\u\\ Ll( n_ [X j {]) . 

By taking limit r — > d~(X) in the above inequality, we derive ( 15.361 . ■ 
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Now we are ready to prove (15.341 . Take R = (R max A \X - Y\&>)/4. We may assume 
that d~(X) < R/2 because otherwise (15.341 follows from ( 15.331 . We then set u to be the 
columns of G(-, Y) in Lemma l5.35l to obtain 

(5.38) \G(X, Y)\ < Cd-{Xf"R- dll - l -^\\G{; F)lk(P-PWD. * = A |X — Y\&)/4. 
Next, we consider the following three possible cases. 

Case 1: \x — y\ < y/t - s < R max . In this case R = y/t - s/4 = \X- Y]gt/4 and thus, we get 
from (EH) and (l3l4t that 

\G(X, Y)\ < Nd-(X)«>R- d/2 - l -«>\\G(; Y)\\i*(P-Ixjcd < Nd-(X)^R- d -^, 

which immediately implies (15.34b in this case. 

Case 2: y/t - s < \x — y\ A R max . In this case R — (\x — y\ A R ma x)/4. We denote Z = (r, z) 
and claim that for all Z 6 Q(X, 2R), we have 

(5.39) \G(r,z, s,y)\ < NC x {t - s y d/2 exp{-K\x -y\ 2 /4(t- s)) , 

where C\ and k are the same constants as in ( 15.331 ) and N = N(d, k). To prove the claim, 
first note that we may assume Y = without loss of generality. Then by ( 15.331 ) we have 

\G(r,z, s,y)\ < C ir - dll e-^^ r X(Q ^{ r ) < C^e^^ Xi^if), 

where we used \z\ — \z — y\>\x — y\/2 = \x\/2. Let us denote 

8 (t) = T- d/2 e-^ 2/4T X(Q ^(T), g (T) = T- dl2 e- K ^x^){T). 

Then the claim ( 15.391 ) will follow if we show that there exists a positive number = N(d, k) 
such that g{r) < Ng{t) for all r < t < \x\~, which in turn will follow if we show that 
go(ri) < Ngo(r2) for all r\ < r2 < 1. But the latter assertion is easy to verify by an 
elementary analysis of the function go. 

We have thus proved ( 15.391 ). which combined with ( 15.381 ) yields 

\G(X, Y)\ < CdrpCfBT^it - sY dl2 exp {-k\x - y\ 2 /4(t - s)} . 

Therefore, we also obtain ( 15.341 ) in this case. 



Case 3: R max < yi - s. In this case R = R max l4, and the desired estimate ( 15.341 ) becomes 
(5.40) \G(t,x,s,y)\ < C{d~ {X) I R max r R; } i x exp{-K|* - y\ 2 /4(t - s)}. 

Since t- s> 16R 2 , for all Z = (r, z) e Q_(Z, 2R), we have 

Iz-yl 2 ! / Jx -y\ 2 /2-\z- x\ 2 \ ^ K/4 __/ K\x-y\ : 



(5.41) exp < -k — \ < exp { -k ■ — ■ } < e K <* exp \ - 

\ r - s J { t-s J { 2{t - s) 

Then, from (15.38b , (15.33b . and d5.4U , we obtain ( 15.40b , which implies ( 15.34b in this case. 

We have thus proved that the estimate ( 15.34b holds in all possible cases. Finally, notice 
that Lemma l535l remains valid if _Sf , X, d~(X), Q-[X, R], and R], respectively, are 

replaced by f _Sf, Y, d + {Y), Q.+[Y,R], and &>C1+[Y,RI Therefore, by replicating the above 
argument to G(-,X), utilizing the estimate ( 15.341 ) instead of ( 15.33b . and using the identity 
(13.4b . we obtain 

\G(t, x, s,y)\ < NS 2 (X, yy%(*. Yf[(f - s) A R 2 ma J- d/2 exp [-k\x - y| 2 /16(f - s)) . 
By replacing k by /c/16, we obtain the desired estimate ( 13.17b . The theorem is proved. ■ 
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6. Appendix 

Lemma 6.1. Assume the condition (LH). Then the condition (LB) is satisfied with the 
same R max and No = No(n, m,/xo, Ni). 

Proof. Let u be a weak solution of J??h = in Cl-[X, R] vanishing on ^Q.-[X,R], where 
X € £2 and R e (0,R ma x)- By using the triangle inequality, for all Y e Q-(X,R/2) and 
Ze Q-(Y,R/2), we have 

|h(F)| 2 < 2\u(Y) - u(Z)\ 2 + 2\u(Z)\ 2 < NR 2 »°[u]l o/24io . Q (XR) + N\u(Z)\ 2 . 

Then by taking average over Z e Q-(Y,R/2) and using (LH), we obtain 



2-d\\ f\\2 

'»L l (fX.lX M' 



where N = N(d, m,po, N\). This proves the part i) of the condition (LB). The proof for the 
other part is very similar and is omitted. ■ 

Lemma 6.2. Let Q.be a time-varying H\ (graph) domain in R rf+1 . Assume that u is a weak 
solution of J£u = f in Q.-[Xo,R] vanishing on &Q.-[Xo,R], where f 6 L^iSl-lXo, R]), 
and denote u — Xci-[X ,R] u - Then we have 

(6.3) f \u - (u) X(hR \ 2 < NR 2 f \Du\ 2 + NR 

JQ-(X a ,R) Jn_[X ,R] 

where (u) Xo ,r - j Q ^ R) u and N = N(d, m, v). 

Proof. We modify the proof of l30l Lemma 3], Without loss of generality, we may assume 
Xq = 0. Let f = g(x) be a smooth function defined on R d such that 

0<£<1, supple B(R), £ = 1 on B(R/2), and \D£\<4/R. 

Setting 6~ x = J £(x) dx and denote 

J3(t):=6 f £(x)u(t,x)dx = 6 f ({x)u(t, x) dx, /3 := R 2 f J3(t)dt. 

JB(R) Joi(t)r\B(R) J-R 1 

Since (u)q r minimizes the integral f \u — c\ 2 among c 6 R"\ we obtain 

(6.4) f |h-(h)o,«| 2 < f \u-~P\ 2 <2 f |fi-£(t)| 2 + 2 f \Pit)-Pt 

Jq-(R) Jq-(R) Jo-Hi) Jo-Hi) 

Notice that u e W^' l (Q-(R)) and Du = xn~[R] Du in Q-(R). Therefore, by a variant of 
Poincare's inequality, we have 

(6.5) f \u -/3{t)\ 2 dX = f f \u{t,x)- P(t)\ 2 dxdt <N f \Du\ 2 dX. 

Jq-(R) J-r 2 Jb(R) Jn_[«] 

We claim that for all s and t satisfying -R 2 < s < t < 0, we have 



(6.6) W)-P(s)\ 2 < NR'" f |Dh| 2 + MT^I/II* M) . 

Jn_[R] 

Assume the estimate ( 16.6b for the moment. By the definition of fi, we then obtain 

f ^8(0-)8| 2 ^ = |B(i?)| f 

JQ-(R) J-r 2 

-0 W) 



<^- 2 f f \B(t)-J3(s)\ 2 dsdt<NR 2 f \Du\ 2 dX + NR 2 ' 11 

J-R 1 J-R 2 Jn_[R] 



L,(n_[R])- 

[«] 
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By combining d6.41 i, (16.51 . and (16.7b . we obtain ( 16.3b . 

It remains us to prove the estimate ( 16.61 ). Setting 17 = i](x) = ij(x)A, where A e R m is a 
constant column vector, and following the calculation in Brown et al. (3), we obtain 



I £(x)u(t, x) ■ A dx - I £ (jc)h(5, x) • A dx 

Ja)(t)nB(R) Joj(s)DB(R) 

+ \ A T A a/1 (X)D B u(X)D a ((x)dX = \ f(X)-A£(x)dX. 



Notice that 6~ l > 2~ d \B{X)\R d . Therefore, by using the properties of the function we get 
for all s and t satisfying -R 2 < s < t < that 

(fi(t) - 0(s)) ■ A < NR~ d ~ l \A\ \ \Du\dX + NR~ d \A\ \ \f\dX. 

By taking A = fl(t) — in the above inequality and using Holder's inequality and 
Cauchy's inequality with e, we obtain (16.61 1. The proof is complete. ■ 

Lemma 6.8. Let Q. be a time-varying H\ (graph) domain in W ,+l . Let a a P satisfy ( 14.5b 
and let S be as in ( 14.6b , where A", are the coefficients of the operator Jzf. Then, there 
exists Sq = £o(d, vq, M) > such that if S < Sq, then the condition (LH) is satisfied with 
po = Po(d,VQ, M), Rmax = R a > and Ni = Ni (d, m, Vq, M). Here, we set R a = 00 if Q. is a 
time-varying Hi graph domain. 

Proof. We shall prove below that there exists a number S\ = S\(d, vo,M) > such that 
if S < S\, then the following holds: There exist positive constants p\ = p\{d,VQ,M) and 
Ci = C\ (d, m, vo, M) such that for any X e 3Q = 0>Q. and R e (0,R a ), if u is a weak 
solution of ^fu = in Cl-[X, R] vanishing on ,^Q.^[X,R], then we have 

(6.9) f \Du\ 2 < Ci (P\ d+2> " f \Du\ 2 , VO < p < r < R. 
Jn-[x,p\ \ r 1 Jn_[x,r] 

We also note that by Lemma 2.2], there is S 2 = S 2 (d, v ) > such that if $ < S 2 , 
then the following holds: There exists a constant p 2 = V-i(d, vo) € (0, 1] such that if u is a 
weak solution of Jzfw = in Q-(X,R) c Q, then we have 

(6.10) f \Du\ 2 < C 2 (P\ d ^ 2 f \ Du \2 f V0 < p < r < R, 

jQ-(X,p) \ r > JQ-(X.r) 

where C 2 = C 2 (d, m, vo). Then we combine (16.9b and ( 16.10b , via a standard method in 
boundary regularity theory to conclude that if S < S\ A S 2 -: Sq, then for all X e D. and 
< R < R a , the following holds: If u is a weak solution of 5£u — in £l-\X, R] vanishing 
on SQ-tX, R], then we have 

(6.11) f \Du\ 2 < N (-)' +2n f \Du\ 2 , V0<p<r<R, 

where po = p\ A /i2 and AT = N(d, m, vq, M). By Lemma [6721 the estimate (16.11b . and the 
energy inequality of Brown et. al (3), we have for all Y e Q-(X, R/4-) and r e (0, R/4-] that 

f |S - (MM 2 < Nr 2 f \Du\ 2 <Nr 2 ( r -] +2 "° f |Dh| 2 

\u\ 2 <Nr d+2+2fi °R- 2 »° f \u\ 2 , 
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where N = N(d, m, vo, M). Then by the Campanato's characterization of Holder continu- 
ous functions (see e.g., JH Lemma 2.5]), we obtain 



JOJX.R) 



Vo/2,jU ;fi-(^/4) 

JQ-(XJt) 

Then, the above inequality together with a standard covering argument yields part i) of the 
condition (LH). The other part of the condition (LH) is similarly obtained. 

Now, it only remains for us to prove the estimate ( 16.91 ). For X e <9Q and R e (0,R„) 
given, let u be a weak solution of Jzfu = in Q_[X, R] vanishing on ^0_[Z,J?]. Denote 
by Jz?o the parabolic operator acting on scalar functions v as follows: 

J? v = v f - D a {a aP D p v). 

For r 6 (0, /?], let v ! be a unique weak solution in Vz(CL.0t, r)) of the problem 

Jz? v'=0 in a_[l,r], 
v'' = m' on 0*(Cl-[X,r\\ 

where i = 1, . . . , m. Existence of such v ! follows from Brown et al. y). We claim that 
there are positive constants = fi(d, vq, M) and AT = N(d, m, vo, M) such that the following 
estimate holds: 

(6.12) f \Dv\ 2 < n(-]' +2>1 f \Dv\ 2 , V0<p<r. 

We may assume that p < r/8 because otherwise ( 16.121 becomes trivial. Since each v ! 
vanishes on r], it follows from (22 Theorem 6.32] and 11221 Theorem 6.30] that 

there exist jj. = fi(d, vq, M) > and N = N(d, v , M) > such that 

(6.13) osc v' <AyV-" sup \vUNffr-^ d l 2 -%\ 7(n ^ r/2]) . 

n -t?,2p] £!_[£, r/4] 

In particular, the estimate (16.13l l implies v'(X) = 0. Then, by the energy inequality of 
Brown et al. [3 1 and [14 Lemma 4.2], we obtain (recall that p < r/8) 



f \Dvf<Np- 2 f |v''| 2 =AVT 2 f |v''(F)-y''(l)| 2 t/y 
Jn_[x,p] Jn_[z,2p] Jn^[z,2p] 

<Np d [ osc v') <N(-) d+2M r- 2 f |v'| 2 

Ul_[X,2p] / \rt Jd.lX.r/2] 



i — 1 , . . . , m. 

1-[X,i 



\ r ' Jn_[l,rl 



where AT = N(d, Vo, M). This completes the proof of the estimate ( 16.121 ). Next, notice that 
w '.— u — v belongs to V2(Q_[X, r]), vanishes on &{Q.-[X, r]), and satisfies weakly 

S&w=D a ({A'*-e!*I m )D0u). 
Therefore, by the energy inequality of Brown et al. |3), we obtain 



\Du\ 2 , 



(6.14) \ \Dw\ 2 <N£ 2 \ 

where <S is defined as in ( 14.61 ). By combining ( 16.121 ) and ( 16.14b . we obtain 

f \Du\ 2 < N (8.\ d+2tl f \Du\ 2 + NS 2 f \Du\ 2 , VO < p < r. 
Jn\x,p] \ r > Ja.\x,r\ Jn-\x,r] 
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Now, choose a number fi\ € (0,/j). Then, by a well known iteration argument (see, e.g., 
ifTTl §111.2]), we find that there exists S\ such that if S < S\, then we have the estimate 
d6.9l >. The lemma is proved. ■ 
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